Paradoxes of magnetorotational instability and their geometrical resolution 
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The magnetorotational instability (MRI) triggers turbulence and enables outward transport of 
angular momentum in hydrodynamically stable accretion discs. By using the WKB approximation 
and methods of singular function theory, we resolve two different paradoxes of MRI that appear 
in the limits of infinite and vanishing magnetic Prandtl number. For the latter case we derive a 
new strict limit of the critical Rossby number. This new limit of Roc = —0.802, which appears for 
a finite Lundquist number of Lu = 0.618, extends the formerly known inductionless Liu limit of 
Roc = -0.828 valid at Lu = 0. 
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The magnetorotational instability (MRI) is the main 
candidate to explain the fast formation of stars and black 
holes by triggering turbulence and angular momentum 
transport in accretion disks. In its standard version 
(SMRI), with a vertical field being applied, the instability 
is non-oscillatory [T-^, while a helical applied magnetic 
field leads to an oscillatory instability (HMRI) [4]. 

Both in the astrophysical context [5] , as well as in lab- 
oratory experiments [6j] , it is vital to know which laws of 
differential rotation are susceptible to MRI. The hydro- 
dynamic reference point is Rayleigh's criterion stating 
that a rotating flow with an outwardly increasing angular 
momentum is stable. This implies, e.g., that a Taylor- 
Couette flow of an inviscid fluid between the inner and 
outer co-axial cylinders of radii Ri < Rq and of infinite 
length that rotate with different angular velocities, Q{Ri) 
and Q{Ro), is stable, if and only if RfQ{Ri) < RlQ{Ro). 
In contrast to this, assuming a perfectly conducting fluid 
and a vertical magnetic field being applied, Velikhov 
[H and Chandrasekhar Q found the more restrictive con- 
dition for stability in the form Q^Ri) < Q^Rq). Remark- 
ably, the latter criterion does not depend on the magnetic 
field strength, i.e. in the limit ^ it does not reduce 
to Rayleigh's criterion valid for B^ = 0. This 'curious be- 
havior of ostensibly changing the Rayleigh criterion dis- 
continuously' [8] constitutes the Velikhov- Chandrasekhar 
paradox which implies a dependence of the instability 
threshold on the sequence of taking the two limits of van- 
ishing magnetic field and vanishing electrical resistivity. 
Its physical origin has been attributed to the fact that 
in a fluid of zero resistivity the magnetic field lines are 
permanently attached to the fluid, independent of the 
strength of the magnetic field [H, [1] . 

Another paradox of MRI emerges in the opposite limit 
of vanishing electrical conductivity. This paradox of in- 
ductionless HMRI [9] refers to the fact that in a helical 
magnetic field a perturbation can grow exponentially al- 
though the instantaneous growth of the energy of any 
perturbation must be smaller than in the field- free case. 

Actually, the astrophysical relevance of HMRI is still 



under debate. On the first glance, according to the 
criterion of Liu et al. it could only work for 

rather steep rotation profiles Q{R) with Rossby numbers 
Ro := R{2n)-^dn/dR < 2 - 2V^ -0.828. This would 
clearly exclude any relevance of HMRI for Keplerian pro- 
files characterized by Ro = —0.75. It has to be noted, 
though, that the validity of the underlying local WKB 
approximation, and the possible role of electrical bound- 
aries for extending the applicability of HMRI to higher 
Rossby numbers are controversially discussed [11]. Sur- 
prising new arguments arose recently from investigations 
of the saturation regime of MRI. For the case of small 
magnetic Prandtl numbers (as they are typical for proto- 
planetary disks), Umurhan speculated about a saturated 
rotation profile with regions of reduced shear, sandwiched 
by regions of strengthened shear [12]. For those latter re- 
gions with steeper than Keplerian profiles, HMRI could 
indeed become of relevance. 

In this Letter we both find the ultimate upper limit 
of the critical Rossby number for HMRI, and resolve the 
mentioned paradoxes. We establish that these physical 
effects sharply correspond to the geometric singularities 
that are inherent on the stability boundaries of leading- 
order WKB equations. 

We start with the local WKB equations for the axisym- 
metric perturbation of a steady-state rotational flow of 
a viscous and resistive fluid in the presence of an axial 
magnetic field that were derived and discussed by several 
authors [H, EH, . They can be rewritten in the typical 
form of a non-conservative gyroscopic system |15l] . 

+ (L> + Qo{l + a^)J)u + (TV + K)u = 0, (1) 

where u = (uR^u^/y)^ is the fluid velocity in polar co- 
ordinates (R^cj)). Separating the time-dependence ac- 
cording to u = uex.-p{jt) yields the eigenvalue problem 
L{j)u = for the growth rate of the perturbation 7, 
where L(7) = 7^/ + 7(L) + (^o(l + c^^) J) + + i^, / the 
2x2 unit matrix, A^ = Qo{ujrj{l + a^) + Ro(cj^ — uJv))J^ 
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with ki2 = fto{un{l — a'^) + Ro(cj^ — o;^)), and 
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In the above equations, = uk'^ and cj^y = r]k'^ are the 
viscous and resistive frequencies, uja = kzB'^{fiop)~^^'^ 
the Alfven frequency, ku and kz the radial and axial wave 



numbers of the perturbation, k 



a 



kz/kj 



= ^{Ro) and Ro = Ko{Ro)^ where Rq is the radial co- 
ordinate of a fiducial point for the local stability analysis. 
We use the convention that p = const is the density of the 
fluid, u = const the kinematic viscosity, 77 = {jiioo-)~^ the 
magnetic diffusivity, a the conductivity of the fluid, and 
po the magnetic permeability of free space. For a = 1, 
cjzy = 0, and cj^^ = 0, Eq. ([T]) is similar to the Hill equa- 
tion for two orbiting mass points connected by a spring 
[Tel, a paradigmatic model of SMRI 0, 

Stable perturbations have 5R7 < provided that 7 
with 5R7 = is a semi-simple eigenvalue of L{j). The 
growing solutions of SMRI are non-oscillatory with ^7 = 
0. Therefore, 7 = implies that det{N + = at the 
threshold of SMRI [3] which results in the critical Rossby 
number (above which the flow is stable) 
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where Re = aftoUJ~^ is the Reynolds number, Pm = 
uJu^~^ the magnetic Prandtl number, and Lu = uja^^^ 
the Lundquist number. Formula dU coincides with that 
derived in [l^ from the Routh-Hurwitz criterion |l7j . 

The Velikhov-Chandrasekhar paradox occurs at infl- 
nite Pm and means that in the ideal MHD case (cj^ = 0, 
cjjy = 0) the limit uoa yields Velikhov's value Roc = 
as the instability threshold rather than Rayleigh's limit 
Roc = — 1 of the non-magnetic case {uoa = 0, cj^y = 0). 

With UOA = s cos ip and uon = e sin (/? in (|4]), we obtain 
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which for £ ^ reduces to 
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Introducing the new parameter Ro' = (1 + 4Re^(l + 
2Ro))(l + 4Re^)-^ we flnd that in the {ujA.^r^.^o')- 
space Eq. ([6j) deflnes a so-called ruled surface {£^(p) ^ 
cos (/:),£: sin cos n(/:)) with n = 2, which is a canonical 
equation for the Pliicker conoid of degree n = 2 [18]. The 
surface according to Eq. (j4]) tends to the Pliicker conoid 

when e = \ oj\ + uj'^ goes to zero. This surface is shown 
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FIG. 1: (a) The critical Rossby number of SMRI as a func- 
tion of UOA ^ LuPm~^ and uor^ ^ Pm~^ for cjjy = 1, a = 1, 
Qo = 1, i.e. for Re = 1. (b) Top view of the surface, (c) 
Cross-sections of the surface along the rays specified by the 
Lundquist number, or, equivalent ly, by the angle Lp that varies 
from to 1.5 through the equal intervals /\(p = 0.1; the hori- 
zontal line corresponds to (p — 7r/2. Note that negative values 
of cjry and e are not physical. 



in the (a;^, ^77, Ro)-space in Fig. [TJa) and in projection 
to the (cja, co'r7)-plane in Fig. [Hb) for Re = 1. For each 
a, ujj^^ and VLq it has the same Pliicker conoid singular- 
ity, i.e. an interval of self-intersection along the Ro-axis 
and two Whitney umbrella singular points at its ends. 
This singular structure implies non-uniqueness for the 
critical Rossby number when simultaneously uja = and 
cj^ = 0. Indeed, for a given Lu, tending the magnetic 



fleld to zero along a ray uja 



^Lu in the {uja^^ti)- 



plane results in a value of the Rossby number specifled 
by Eq. (j6j), see Fig. [U^c). The limit value of the critical 
Rossby number oscillates between the ideal MHD value 
Roc = for Lu = 00 ((/:) = 0) and the non-magnetic 
value Roc = — 1 — (2Re)~^ for Lu = {^p = tt/2), which 
explains the Velikhov-Chandrasekhar paradox. 

Now we turn over to the paradox of induction- 
less HMRI which is related to a similar geometric 
singularity as discussed above. The leading order 
WKB equations that describe the onset of instabil- 
ity of a hydrodynamically stable TC-flow with a he- 
lical external magnetic fleld are ^ = with = 
{ur, u^, Br{iiqp)-^/'^,B^{iiqp)-^/'^) and 
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where the additional parameter uja^ = Rq^ B'^{pop)~^^'^ 
is the Alfven frequency of the azimuthal magnetic fleld 
component [14]. For uja^ = these equations yield ([1]). 
The dispersion equation det(i^ — 7/) = reads 

A"^ + aiX^ + a2A^ + (as + ibs)X + a4 + ib4 = 0, (8) 
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where / is a 4 x 4 unit matrix, A = 7(cjj^a;^) and 



2(1 + Pm"^)yPm, 



(9) 



ai 

a2 = 2(l + (l + 2/3^)Ha^)+4Re^(l + Ro)Pm + a^/4, 

as = ai(l + (1 + 2/3^)Ha^) +8Re^(l + Ro)VPm, 

a4 = (l+Ha^)%4/32Ha^+4Re^(l+Ro(PmHa^ + l)), 

63 = -8/3Ha^ReV^, 64=fe(l+(l-Pm)Ro/2)/VPm, 

where we have introduced now the Hartmann num- 
ber Ha = LuPm~^/^ and the hehcity parameter j3 = 
aujA^^^ of the external magnetic field. According to the 
analogue of the Routh-Hurwitz conditions for the com- 
plex polynomials — the Bilharz criterion [l7] — the thresh- 
old of HMRI is defined by m4(/3, Re, Ha, Pm, Ro) = 0, 
where 7714 is the determinant of the so-called Bilharz ma- 
trix [l^, [13, [HI composed of the coefficients (|9]). The 
stability condition 5RA < holds if and only if 7714 > 
[13, [l^^ [HI • For /3 = the dispersion equation and thus 
the threshold for HMRI reduce to that of SMRI [li]. 

In the following we will see that in the limit Pm 
it is again Lu that governs the value of Roc- For this 
purpose, we show in Fig. [2fa) a typical critical surface 
m4 = in the (Pm,Re"\Ro) -space for the special pa- 
rameter choice Ha = 15 and (3 = 0.7. On the Ro-axis we 
find a self-intersection and two Whitney umbrella sin- 
gularities at its ends. At the upper singular point, i.e. 
exactly at Pm = 0, we get (see [lj|) 
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In the limit Ha — )• oo, this critical value is majorated by 



Roc(/3) 



1 + 4/?^- (1 + 2/32)^1 + 4/34 
2/?2 
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with the maximum at the well-known Liu limit Roc = 
2 - 2^2 - -0.828 when /3 = V2/2 0.707 [Toj, [Hi. 

In Fig. [2fa) we see that the case with Pm = is con- 
nected to the case Pm ^ by the Pliicker conoid singu- 
larity, quite similar as it was discussed for the paradox 
of Velikhov and Chandrasekhar. Interestingly, Roc for 
the onset of HMRI can increase when Pm departs from 
zero which happens along curved pockets of HMRI, see 
Fig.[2fa). The two side bumps of the curve Re-^(Pm) m 
a horizontal slice of the surface correspond to the domains 
of the essential HMRI while the central hill marks the he- 
lically modified SM/?/ domain, according to the classifica- 
tion introduced in [1^] . For small Pm the essential HMRI 
occurs at higher Ro than the helically modified SMRI, 
while for some finite value of Pm the central hill and the 
side bumps get the same value of Roc . Most remarkably. 
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FIG. 2: (a) The critical Rossby number of the essential HMRI 
and helically modified SMRI for Ha = 15 and P = 0.7 (left) 
m the (Pm,Re"\Ro) -space, (b) The critical Rossby num- 
ber for Lu = 0.5 and /S — 0.6 in the (Pm, Re~^, Ro)-space 
and (c) its cross-sections in the (Pm, Re~^)-plane for (black) 
Ro=-0.842, (blue) Ro= -0.832, (green) Ro=-0.822, (red) Ro=- 
0.812, (brown) Ro=-0.802. 



there is a value of Roc at which the two side bumps of 
the curve Re~'^(Pm) disappear completely. This is the 
maximal possible value for the essential HMRI, at least 
at the given p and Ha. Now we can ask: how does this 
limit behave if we send Ha to infinity, and to which value 
of Lu does this correspond? 

Actually, with the increase in Ha the stability bound- 
ary preserves its shape and simultaneously it com- 
presses in the direction of zero Pm. Substituting Ha = 
LuPm~^/^ into the equations ([9]), we plot again the sur- 
face 7714 = in the (Pm, Re~^, Ro)-space, but now for a 
given P and Lu, Fig. [Sfb). 

The corresponding cross-sections of the instability do- 
main in the (Re~^, Pm)-plane are shown in Fig.[2fc). At 
a given value of Ro there exist three domains of instabil- 
ity with the boundaries shown in blue and green. Two 
sub-domains that have a form of a petal correspond to 
the HMRI. They are bounded by closed curves with a 
self-intersection singularity at the origin. They are also 
elongated in a preferred direction that in the (Re~^, Pm)- 
plane corresponds to a limited range of the magnetic 
Reynolds number Rm = PmRe. The central domain, 
which corresponds to the helically modified SMRI, has a 
similar singularity at the origin and is unbounded in the 
positive Pm-direction. In comparison with the central 
domain, the side petals have lower values of Rm. 

Now we reconsider again the limit Pm 0, while 
keeping Lu as a free parameter. At the origin all the 
boundaries of the petals can be approximated by the 
straight lines Pm = RmRe~^. Substituting this ex- 
pression into equation 7714 = 0, we find that the only 
term that does not depend on Pm is a polynomial 
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FIG. 3: (a) Discriminant surface in the (Lu, /3, Ro)-space and 
(b) its cross-section at /3 = 0.634. (c) Interaction parameter 
N = Lu^Rm"^ at the essential HMRI maxima. 



(5(Rm, Lu, /3, Ro) = po-\-piRm^ -\-p2Rm!^ -\-psRm.^ , where, 
e.g., po = Lu'^(4/34Lu^ + 2/32 + 4Lu^/32 + l)^ 

The roots of the polynoraial are coefficients Km of the 
hnear approximation to the instabihty domains at the 
origin in the (Re~^, Pm)-plane. Simple roots mean non- 
degenerate self-intersection of the stability boundary at 
the origin. Double roots correspond to a degeneration 
of the angle of the self-intersection when it collapses to 
zero which happens only at the maximal critical Rossby 
number, Fig.lSfb). In the (Lu, /3, Ro)-space a set of points 
that correspond to multiple roots of the polynomial Q is 
given by the discriminant surface 64A^poP3 = [19]. The 
surface ps = consists of a sheet Ro = — (1 + Lu^)~^ 
corresponding to the doubly degenerate infinite values 
of Rm at the maxima of the helically modified SMRI. It 
smoothly touches along the /3-axis the surface A = that 
consists of two smooth sheets that touch each other along 
a spatial curve — the cuspidal edge — corresponding to 
triple roots of the polynomial Q, Fig. [H^a). 

Every point on the upper sheet of the surface A = 
represents a degenerate linear approximation to the es- 
sential HMRI domain and therefore a maximal Ro at the 
corresponding values of /3 and Lu. Numerical optimiza- 
tion results in the new ultimate limit for HMRI Roc — 
-0.802 at Lu ^ 0.618, (3 0.634, and Rm 0.770, see 
Fig. [3l^b). This new limit of Roc on the cuspidal edge is 
smoothly connected to the inductionless Liu limit by the 
upper sheet of the discriminant surface, which converges 
to the curve (pTj) when Lu = 0. We point out that the 
new limit is achieved at Ha 00 when the optimal Pm 
tends to zero in such a way that Lu 0.618. Figure [31(c) 
shows the behaviour of the so-called interaction param- 
eter (or Elsasser number) N = Lu^/Rm for the HMRI 
sheet. It is remarkable that, at Lu = 0, HMRI starts to 
work already at N = 0. This can be explained by the ob- 
servation that the optimal value for HMRI corresponds 
to NHa = LuV(RmV^) = 1/(1 + 2'^/^) = 0.586, [M]- 



Later, for increasing Lu, the optimal N acquires final val- 
ues, passes through its maximum and at Lu 0.618 and 
/3 0.634 it terminates at N = 0.496. 

Inspired by the theory of dissipation induced instabili- 
ties [15| , we have resolved the two paradoxes of SMRI and 
HMRI in the limits of infinite and zero magnetic Prandtl 
number, respectively, by establishing their sharp corre- 
spondence to singularities on the instability thresholds. 
In either case, it is the local Pliicker conoid structure that 
explains the non-uniqueness of the critical Rossby num- 
ber, and its crucial dependence on the Lundquist number. 
For HMRI, we have found an extension of the former Liu 
limit Roc — —0.828 (valid for Lu = 0) to a somewhat 
higher value Ro —0.802 at Lu = 0.618 which is, how- 
ever, still below the Kepler value. To study the possible 
consequences of this new limit for the saturation of MRI 
in accretion disks or experiments, is left for future work. 
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